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We consider the effect of broadband decoherence on the performance of refocusing sequences, 
having in mind applications of dynamical decoupling in concatenation with quantum error correcting 
codes as the first stage of coherence protection. Specifically, we construct cumulant expansions of 
effective decoherence operators for a qubit driven by a pulse of a generic symmetric shape, and 
for several sequences of n- and 7r/2-pulses. While, in general, the performance of soft pulses in 
decoupling sequences in the presence of Markovian decoherence is worse than that of the ideal 
5-pulses, it can be substantially improved by shaping. 



I. INTRODUCTION 




Dynamical decoupling |l], ||, ||] (DD) can be very effec- 
tive in protecting coherence of a quantum system against 

. This, 



low-frequency environment! 

H g g§ gj, 0, g ^ |2c 

combined with low resource requirement, makes it attrac- 
tive as the first-level coherence protection technique, in 
combination with quantum error correcting codes |27| ] 
(QECC). For such a combined decoherence protection 
technique to work universally, the performance of the DD 
should not be reduced in the presence of fast environment 
modes whose effect on coherence is to be dealt with by 
QECC. 

Previously, the effects of broadband noise on DD were 
analyzed in a number of publications, with the primary 
target being the 1// or telegraph noise |l6|, ^8|, [2{j, 
|30| , |3j| . However, the effectiveness of refocusing se- 
quences was mostly studied numerically, apart from spe- 
cial exactly-solvable cases |3(], ^2j. No attempt has been 
made to investigate general properties and limitations of 
decoupling under these conditions. Certainly, there were 
no analytical studies of effects of pulse shaping on decou- 
pling in the presence of broadband decoherence. 

In this work, we concentrate on the dynamical decou- 
pling (DD) of a single qubit (spin) in the presence of 
a broad-spectrum oscillator bath [Fig. We first con- 
sider the effect of the pulse shape in the special case 
where the external spin couplings can be described by 
a combination of a Markovian noise (described by the 
Lindblad decoherence operators) and a time- independent 
magnetic field. This basic problem is similar in impor- 
tance to the canonical pulse-shaping for a nuclear spin 
in the presence of a chemical shift. For a qubit driven 
by an arbitrarily-shaped pulse, we construct an analyt- 
ical expansion of the average decoherence operator, an 
analogue of the average Hamiltonian expansion Q 
but for the qubit density matrix evolution. Analyzing 
the first two terms of such expansions for several decou- 
pling sequences, we formulate the conditions necessary 



for improved coherence of the qubit, and construct nu- 
merically the pulse shapes that satisfy these conditions. 
We then compare the performance of the obtained pulse 
shapes with "hard" pulses and the conventional first- and 
second-order nmr-style self-refocusing pulses in several 
decoupling sequences both analytically and numerically. 
For numerics, we model the oscillator bath as a com- 
bination of classical correlated gaussian noise and the 
Markovian noise [Fig. 0(a) and |l|(c) respectively] . 
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FIG. I: Assumed spectral function of a generic oscillator bath 
(schematic). Dynamical decoupling can be very effective in 
preserving coherence in the presence of (a) featureless low- 
frequency bath[^f[ and also (b) sharp resonances J35J]. In this 
work we consider the effect of Markovian noise (c) on DD of a 
single qubit. We assume F(uj) =const up to a large cut-off fre- 
quency w max = A, so that the qubit dynamics associated with 
the modes (c) can be modeled by the Linblad equation (p"f|). 
When the high-frequency spectrum is not flat but is, e.g., a 
power law oc \/uo a as in Refs. ^o|, [j^, the decomposition into 
low-frequency (a) and Markovian (c) parts becomes approx- 
imate, with the Markovian decoherence rates 7 determined 
by F(ui) near u> = SI = 2n/r determined by the sequence 
duration r. 

We analyze two kinds of decoupling sequences, the 
usual decoupling sequences utilizing 7r-pulses, and the 
sequences using composite pulses constructed of three 
7r/2 pulses similar to those used in the NMR WaHuHa 
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experiment Q]. While the former sequences interact rel- 
atively little with the decoherence operators (<5-pulses 
leave the diagonal parts of the decoherence operators in- 
variant), the latter ones are constructed to symmetrise 
the decoherence operator between all three channels. 
Such a redistribution seeks to minimize the detrimen- 
tal effects of the non-Hermitian evolution caused by the 
decoherence operators on higher-level control sequences. 



A schematic plot of a single representative component of 
this function is shown in Fig. [[J. 

For simplicity, we assume a combination of weak bath 
coupling and finite temperature T. so that the (generally 
non-Mar kovian) master equation (2l| |3(| [3?], g8[ |3|] in 
the leading order ( "Born approximation" ) is satisfied for 
the combination of Hs + Hcit). For the case of a single 
qubit, which is primarily discussed in this work, 



II. MODEL 



Hamiltonian 



p = -i[H s + H c (t),p}-^J dt'Tl 
-- A fdt'Tf(t-t')[<j a ,{a^p}], 



I T af3 f , _ ,l\ra \R 



(8) 



A system of individually-controlled qubits in the pres- 
ence of an oscillator thermal bath can be described by 
the following idealized Hamiltonian 



H = H c (t) + H s + H SB + H L 
where the control Hamiltonian 



(1) 



(2) 



the system Hamiltonian with one- and few-qubit terms, 



the linear oscillator bath, 

and the bath-coupling Hamiltonians, 



h sb = \ E ^ CT » + 7 



E-a/3 a 

n,n' ,a,/5 



(4) 



(5) 



Here a", a = x,y,z, are the Pauli matrices for n-th 
qubit, a M , are the annihilation and creation operators 

for the oscillators of the bath, and the operators 6", 
describe the various couplings with the bath. A linear 
in phonons coupling, e.g., for the single-spin term can be 
written as 



fa a + p a t 
(2m^) 1 /2 



(0) 



In the limit where the oscillator modes have a continu- 
ous spectrum, the kinetics of the decoherence effects is 
determined by the set of spectral coupling functions (ma- 
trices), e.g., for the single-qubit coupling, 

^ jot 

F n ~ n >a{u) = 77 y~] — n ^ 5{uj - wA (7) 



where the dissipation kernel T\ (t) + iJF(t) = f(t) is de- 
fined as 



F(t)= I -[%)K + l)e' ut + F*(, 



ZUJt ' ^* f uj)n u e- iu> % (9) 



n u = exp((3u>) — 1 is the oscillator occupation number, 
and (3 = ti/T. 

Dynamical decoupling in the presence of a low- 
frequency bath [Fig. 0(a), (b)] was analayzed in Refs. 
|35| . General conclusion is that a carefully- designed se- 
quence with the period r = 27r/f2 provides an excellent 
decoherence protection as long as the adiabaticity con- 
dition £1 > lo c is satisfied. Here uj c is the bath cut-off 
frequency such that all lo^ < oj c in Eq. (0). 

More precisely, the analysis in Ref. ^ was done for 
a featureless low-frequency bath [Fig Ua)] in the ap- 
proximation of non-Markovian quantum kinetic (master) 
equation. The approach was to design a decoupling se- 
quence effective for a closed system with the "frozen" 
bath, where each term in the bath-coupling Hamiltonian 
Hsb [Eq. (^)] is replaced by the corresponding term with 
a non-zero c- number coefficient in Eq. (|^) [e.g., 6" — > B", 
Jn'n' ~~ ¥ Jn'n'' e ^ c -' with the matching indices.] The cor- 
responding unitary evolution operator U(t — r) after one 
decoupling period can be characterized by the decoupling 
order K, the number of terms in the Magnus (cumulant) 
expansion in powers of Hs which are suppressed identi- 
cally. Such a suppression can be expressed in terms of 
certain algebraic conditions which were used in the anal- 
ysis of the solutions of the full master equation, order- 
by-order in powers of the adiabaticity parameter ui c /Q. 
At K — 1, the single-phonon T\ processes are completely 
suppressed as long as w c < O (this can also be seen b y a n- 
alyzing the absorption spectra of the driven system poj). 
Additionally, since the bath coupling is modulated at the 
frequency fl, the effective dephasing rate is reduced by 
the factor of o; e /Opl]|. At K = 2, in many cases (e.g., 
for a single qubit), dephasing rate is suppressed faster 
than any power of u> c /ft < 1: terms of every order in the 
corresponding expansion are suppressed pi|. 

The description in terms of the master equation fails 
in the presence of sharp resonances [Fig. [3(b)]; the cor- 
responding modes have to be included in the Hamil- 
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tonian of the system and the controlled dynamics re- 
analyzed^]. Compared to qubits-only systems, presence 
of oscillators in H$ can generate some additional terms in 
the Magnus expansion of the unitary evolution operator 
U(t) for the closed system which requires more careful 
sequence design. When this is done, only the effective, 
DD-renormalized coupling to such modes matters. Non- 
equilibrium effects like mode heating do not occur as long 
as this renormalized coupling is small compared to either 
the width of the spectral peak or the corresponding fre- 
quency biasp5[. 



B. Markovian decoherence 

We now consider the dynamical decoupling in the pres- 
ence of broadband bath modes, see Fig. [l](c), assuming 
that the corresponding inverse correlation time A = w max 
is much larger than the control bandwidth which deter- 
mines the pulse rate which can be implemented for the 
system. Out of the master equation (ph, we separate a 
part Markovian on the time scale slow compared 

to k/ max , and assume that the remaining relatively slow 
modes are limited to frequencies lo < uj c <C w mm . Most 
generally, Markovian evolution implies the dissipator in 



the Lindblad form 1 41 



A 



(10) 



where p is the density matrix and Aj are the decoherence 
operators. Such a form can also be recovered from Eq. (j^) 
in the Markovian limit where the time-dependence of the 
kernel F{t) becomes (5-function-like. We note that this 
separation of the approximately Markovian part of the 
bath, Eq. (fL0|), does not imply the consitions of high 
temperature or weak coupling combined with adiabatic- 
ity which are usually necessary for the applicability of 
the Markovian master equation |57| |9[ ^2|. 

For example, with flat part of « const as in 

Fig. [|(c), we could set the low-frequency cut-off ui c > 
and include the effect of (2n u + l)fi(u>) in Eq. (^) in the 
non-Markovian part of the bath kernel. This is precisely 
the case solved numerically in Sec. IV where the slow 



bath modes are replaced by correlated classical noise as 
appropriate for fku < T . 

To start, we will follow the approach of Ref. |2l], and 
begin by analyzing the dynamics of the system with 
"frozen" dynamics of the slow modes in Fig. 0(a). That 
is, we will assume that the decoupling frequency Q is 
large compared to u> c , so that any bath couplings can be 
approximated by appropriate terms in H g. Then, the 
effect of the fast environmental modes [Fig. 0(c)] on dy- 
namical decoupling can be analyzed by considering the 
driven Lindblad equation, 



For a single qubit, parametrize the dissipator in terms 
of a Hermitian non-negative-definite rate matrix j a p, 

V(J>) = \ E %p([^ a P,^} + [a a ,p^}). (12) 

It is convenient to separate the symmetric real and anti- 
symmetric imaginary parts of the rate matrix, 



la/3 = lap + iCapylZj. 



(13) 



Then, if we rewrite the density matrix in terms of the 
Bloch vector R = R(i), 

p(t) = \ (1 + R(t) • 8) , R 2 <1, (14) 

the equation for spin kinetics in the presence of the single- 
spin control ^ and system (|J) Hamiltonians becomes 

R = [V(i)xR]-fR-K, (15) 
-f R = [B x R] + (7- lTr 7 )R(i). (16) 

The terms with V(t) and B describe precession in the 
effective magnetic fields of the control and the system 
Hamiltonians, the term with 7 describes the coherence 
loss due to Markovian bath, while the last term in 
Eq. (|l5| ) ensures the correct equilibrium value of R in 
the absence of control. 

If we work in the basis associated with the dissipator, 
the rate matrix 7 is diagonal. In particular, for the case 
of nuclear magnetic resonance (or any physical qubit with 
large energy difference between the levels) , we have j xx — 
= 7, j zz = 70, and the matrix ( |l6| ) takes the form 



7 + 70 B z 



r = 



-B,, 



-B z 7 + 70 B x 



(17) 



B„ 



Brr 



27 



-i[H c (t) + H s ,p]+V{p). 



(11) 



These decay rates correspond to the usual coherence 
times T^ 1 = 27, T^ 1 — 7 + 7^. The corresponding 
vector 1Z has only one component, 1Z Z > 0. In the spe- 
cial case 70 = 7, the rate matrix is proportional to the 
identity matrix, 7 = 7!. 

We emphasize again that, while the analytical analysis 
is done assuming time-independent "magnetic field" B 
in Eq. ([H]), the decoupling accuracy for obtained pulse 
shapes is also verified numerically in Sec. IV with the 
components B^t) taken as a classical time-dependent 
Gaussian correlated noise. 



C. Interaction representation & Magnus expansion 

In dynamical decoupling, it is the control Hamiltonian 
Hc(t) that dominates the dynamics. The evolution as- 
sociated with the system Hamiltonian H$ [Eq. (P)] or 
the dissipator V(p) [Eq. ([To|)] is much slower. Thus, it 
is convenient to consider the dynamics in the interaction 
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representation with respect to Hcif)- To this end, we 
introduce the non-perturbed unitary evolution operator 
U Q = U (t), 



U (t) 



(18) 



where T t represents the usual time-ordering operator, 
and the corresponding orthogonal spin-rotation matrix, 
Qo = Qo(t), can be computed as follows 



U Q o"Ul = Qfo*. 



(19) 



The components Qq satisfy the Bloch equation 

Qf(t)=e^ s V^t)Q 5 f(t), Q (0) = 1. (20) 

If we write the solution of the uniform version of 
Eq. (|l5| ) (i.e., with TZ = 0) in terms of the full evolu- 
tion matrix Q = Qit), 



R(i) = Q(t) R , 



(21) 



we can introduce the decomposition Q = Qo S, where the 
matrix S defines the slow evolution of the Bloch vector 
in the rotating frame defined by the control fields. The 
equation for S reads 



S = -t(t)S, f(t) = Ql(t)tQ (t), 



(22) 



where f (t) is the evolution operator fll6|) in the interac- 
tion representation. The formal solution of Eq. ( p2| ) can 
be again written as a time-ordered exponent, 



S(t) = 



T t exp(-^ di'f(i'))- 



(23) 



For a periodic control field, H c (t + r) = H c (t), such 
that the zeroth-order rotation matrix is also periodic, 
Qo(t + t) = Qo(t), the time-ordered exponent can be 
evaluated in terms of the average decoherence operator 
r, an analogue of the average Hamiltonian [p3|, B4|: 



S(ut) 



-utT 



r = r 



(0) 



-(i) 



where 



r (o) = 1 



(24) 

(25) 
(26) 



(fc) 



etc. Generally, the term T , k > 0, of this expansion 
contains a (k + l)-fold integration of commutators of the 
rotating-frame decoherence operator T(tj); it scales as 

oc r k . Trace of a commutator is zero, thus TrT 
for all k > 0; this implies 



(k) 







TrT = TrT (0) = Trf = 2Tr7. 



(27) 



Generally, the trace-conserving modifications of the ef- 
fective decoherence operator correspond to (some degree 
of) symmetrization between the decoherence channels. 
More precisely, the real part of the smallest eigenvalue of 
r cannot be smaller than that of T. This can be seen by 
separating out the trivial part of T proportional to iden- 
tity matrix, which corresponds to uniform decoherence 
for all channels and remains unchanged in the interac- 
tion representation (22). 



D. Coherence loss in controlled system 



As would be expected on general grounds, Eq. ( |27j ) im- 
plies that the DD cannot eliminate the Markovian dissi- 
pation; at best we can hope to eliminate the off-diagonal 
terms of matrix T and redistribute the rates over the 
decoherence channels. For example, for an NMR qubit 
decoherence model, Eq. (^), redistribution (symmetriza- 
tion) over all three channels would lead to 



r = §<2 7 . 



70)1- 



(28) 



With 7 <S 70, this is equivalent to a 33% reduction of 
the maximum decoherence rate, or a 50% increase of the 
coherence time measured by the fidelity minimized over 
the initial conditions. 

We should note that the quoted estimate for the max- 
imum decoherence rate improvement does not take into 
account the antisymmetric part of the rate matrix 7 
[Eq. (|13|) 1, or, equivalently, the non-uniform term TZ in 
Eq. (p_5|). In the absence of control fields, this term is re- 
sponsible for asymptotic thermal-equilibrium orientation 
of the spin; it is small and can be ignored at sufficiently 
high temperatures (e.g., room-temperature NMR), but 
it can be large at temperatures small compared to the 
qubit level difference A_E i = |B|. 

We can drop the non-uniform term TZ in the analysis 
of decoherence since it does not affect the evolution of 
the average fidelity, that is, the fidelity averaged by the 
initial conditions, T(t) = Tr(p(0)p{t)) o . Indeed, Eq. (|l|) 
is a set of linear differential equations for the components 
of the Bloch vector R(£); the solution with given initial 
conditions R(0) = Ro can be written as a sum of the so- 
lution (21) of the uniform equation with the same initial 
condition, and that of the non-uniform equation but with 
zero initial condition. In the expression for the average 
fidelity, 



^(f) = ^(i + (R -R(t)>R ; 



(29) 



the solution of the uniform equation is bilinear in the 
components of Ro while the non- uniform part is linear; 
the averaging over all directions of R leaves only the bi- 
linear part, 



T(nr) = \ + 1 Tr Q(nr) = \ + \ Tr 



(30) 
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where we assumed the zeroth-order rotation matrix to be 
periodic, Qo{nr) = i, and used the representation j2^ ) 
in terms of the average decoherence operator T. 

The obtained expression (|30|), along with the proper- 
ties of the average decoherence operator T discussed in 
the previous section [see Eq. (E7|)], imply that the aver- 
age fidelity cannot be improved by means of dynamical 
decoupling beyond eliminating the effect of the system 
Hamiltonian H s [Eq. (g)]. With B = 0, this follows di- 
rectly from the convexity of exponentials in Eq. (|3C| ) if 
one works in the basis where the average decoherence 
matrix 7 = diag(7i, 72, 73), < 71 < 72 < 73- In par- 
ticular, the maximum possible fidelity with given Tr7 
corresponds to only one non-zero component, 73 = Tvy. 
In this case two (or any even number) qubit errors com- 
pensate each other. On the other hand, the symmetric 
case 71 = 72 = 73 = TY7/3 in the absence of error cor- 
rection corresponds to the fidelity minimum. 

We should also note that a symmetrization of deco- 
herence rates over channels [complete as in Eq. (^8|), or 
partial] would prohibit the use of special QECCs de- 
signed for strongly asymmetric error rates between the 
channels g§ (E§ g g |47|, ||. 

In this work we assume that the benefits of dynamical 
decoupling (due to reduced decoherence coming from the 
slow degrees of freedom) outweigh the detrimental effects 
related to symmetrization of the coherence rates. 



III. PULSE AND SEQUENCE 
CHARACTERIZATION 

A. Hard pulses 

Let us now use the obtained formalism to analyze con- 
trolled dynamics. We start with hard pulses, with the 
5-function centered at the middle of the interval of du- 
ration t p . For example, for a (4>o) x pulse (rotation angle 
4>q with respect to a;- axis), the non-zero control field in 
Eq. (g) can be written as 

V x (t) = cf>oS(t-T p /2), 0<t<r p . (31) 

For the special case (f> = tt (pulse ir x ), the corresponding 
rotation matrix (EOh is 



X = ir x pulses (which is the repeated part of the Carr- 
Purccll sequence), the first-order average decoherence op- 



Qo(t) 



1, < t < Tp/2; 

diag(l, -1, -1), T p /2<t<T p . 



(32) 



Over the second half of the interval, the components T xy 
and T xz are inverted; as a result in the leading order 



lyy + Izz 



The expression for r 1 (tt x ) is too complicated to emote 
here. However, for the symmetric sequence 2s = XX of 



Izz + Ixx B x ~i yz ] . (33) 

~B X — Jy Z "f XX + Jyy 



crator f (1) (2s) 



(2s) 



0, while the zeroth-order term V 
is the same as for a single pulse given by Eq. (|33|). 

Deriving similar expressions for a tt v pulse, we 
can easily construct the expressions for standard two- 
dimensional decoupling sequences of hard 7r-pulses. In 
particular, the group-averaging|l^ 4-pulse sequence 
4p = XYXY results in cancellation of all off-diagonal 
terms to leading order, 



r (0) (4 P) 



lyy + Izz 

l xx +l zz I , (34) 

J xx + lyy 



while the symmetric sequence 8s = XYXYYXYX (sym- 
metrized version of 4p) results in cancellation of all off- 
diagonal terms [r (0) (8s) = f (0) (4p), see Eq. @] and 

does not produce any lst-order corrections, r^^s) = 0. 

We note that with the hard pulses, the contribution to 
(2) 

r (8s) is non-zero already for a Hamiltonian evolution 
with 7 = 0, as long as both B x and B y are non-zero. 

The hard tt pulses do not modify the structure of the 
diagonal part of the evolution matrix T. A redistribution 
(symmetrization) over different decoherence channels can 
be done if we deliberately orient the spin along differ- 
ent axes with 7r/2-pulses. In the absence of the system 
Hamiltonian (B = 0), for NMR decoherence model (fT^), 
the symmetrization can be achieved with the WaHuHa 
(MREV-4) NMR experiment Q |50|, a sequence of vr/2 
pulses with the cycle 5 = XY0YX of duration r = 5r p , 
where each pulse is centered as in Eq. (|3l]) and stands 
for a free evolution interval equal to the pulse duration 
Tp. The averaging of the second-rank coupling achieved 
by this sequence is used to eliminate dipolar coupling 
between nuclei. However, the accuracy achieved by this 
sequence is rather sensitive to chemical shifts; the cor- 
rections are present already in the leading order effective 
evolution operator. 

To achieve both the decoherence symmetrization and 
the decoupling, we constructed several pulse sequences 
based upon the composite pulses R a = XYX, Rp 
X Y X . and the inverted versions R a = XYX, and Rp 
X Y X Here, X and Y now denote the 7r/2 pulses applied 
in the corresponding directions. The composite pulses 
R a and Rp are mutually orthogonal 7r-pulses constructed 
in such a way that the spin's x, y, and z axes spend 
equal time oriented along the z-axis. As a result, they 
achieve a leading-order symmetrization of the NMR de- 
cohernce matrix, see Eq. ([l7]) with B = 0. The twelve- 
and twenty-four pulse sequences, 12 and 24 respectively, 
are merely the universal decoupling sequences 4p and 8s 
constructed from these composite pulses, 

12 = R a RpR a Rp, (35) 
24 = R a RpR a Rp RpR a RpR a , (36) 

while the 48-pulse sequence includes an additional cycle 
of phase ramping analogous to that in MLEV-16 and 
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higher order sequences Q, 

48 = R a RpR a Rp R a Rj3R a Rp + (reverse). 



(37) 



The basic 12-pulse sequence achieves symmetrization but 
not the decoupling; the leading-order average decoher- 
ence operator reads 



w(0) 1 



(> 



4(2 7 + 70) 

~2By 

B 7 



2 By 

4(27 + 70) 
B, 



-B z 
-B z 
4(2 7 + 70)- 



(38) 



The longer sequence 24 achieves both the decoupling and 
symmetrization of the decoherence operator (|l7]) in the 

leading order, i.e., T^l is given by Eq. (p8j). The sequence 

48 achieves symmetrization to subleading order, i.e., r 48 

r 0. Moreover, in the 



is given by Eq. ( |28j ) while T 48 
experimentally relevant case B x 

has non-zero elements oc -B, (7— 7^)7^ along the diagonal. 



B y = 0, the r 4 2 8 } only 



With non-zero decoherence matrix 7^0, the evolution 
operator in the interaction representation, Eq. (p2|), will 
contain terms that depend both linearly and bi-linearly 
on the components of the matrix (^). In other words, 
the time-dependence of T(t) will be only through terms 
proportional to cos 0, sin 0, cos 20, and sin 20. The func- 
tions of the doubled angle are specifically due to the 
Markovian decoherence operators; they were not present 
in the analysis of Hamiltonian dynamics in Refs. 35. 5^. 

For the symmetric pulse shape, 

V x (t p -t) = V x (t) , 0(r p - t) = O - 0(i) , (42) 

it is convenient to introduce the symmetrized angle 

<p(t) = <f>(t) - <j> /2, <p(r p -t) = -<p(t). (43) 

The averages of time-dependent terms can be then writ- 
ten as follows: 



B. Finite-duration pulses 

Let us now analyze the performance of discussed se- 
quences when finite-duration "soft" pulses are used in- 
stead of the ideal hard pulses. Basically, we extend the 
formalism of Refs. |35|, |5^ to the case of non-Hcrmitian 
evolution (|l5|) with the most general symmetric rate ma- 
trix 7. 



1. Average evolution operator for a single pulse 

We begin with the qubit evolution driven by a one- 
dimensional pulse of arbitrary shape, 



(COS ( 



9o 1 * i\ .90 
vcos— , (sin 0) = v sm ■ 



(39) 



The (zeroth-order) evolution operator due to control field 
alone is given by 



U c (t) = exp(-KT x 0(t)/2), y(t) = dt'V(t'). (40) 

Jo 

This is just a rotation by angle <fi(t) around the x-axis, 



1 

Q (t) = ( cos <j>(t) sin0(i) 
— sin0(i) cos4>(t) 



(41) 



Specifically for a symmetric ir x pulse, we get 



2 ' x r/ 2 ' 

(cos 20) = v-i cos 4>q , (sin 20) = v 2 sin O > 



where 



v = (cos ip(t)) , v 2 = (cos2(p(t)), 



and the pulse-averages are defined simply as 



1 



(f(t)) = - dtf(t) 



T, 



(44) 
(45) 



(46) 



(47) 



v Jo 



We see that, to leading order, the pulse-shape depen- 
dence of the average evolution operator [Eq. (p5|)1 is 
reduced to only two constants. The parameter v gives 
an effective pulse length to first order; it is equal to 
zero, v — for ideal ^-inversion pulses (rotation angle 
0o = 7r), as well as for Hermitian pulses |33| or other 
lst-order self-refocusing pulses|ll| |5^, pi| . The parame- 
ter v 2 vanishes for ideal <5-pulses with the rotation angle 
0o = tt/2; for hard inversion pulses [Eq. (|3l"|)l we have 
ip = ±7r/2, thus cos 2ip = —1 throughout the interval, 
which gives v 2 = - 1- For soft pulses with limited am- 
plitude, \v2\ < 1. The values of the parameters for some 
other pulse shapes are listed in Table Q. 



r (7r a 



lyy + Izz 
v(-B y +j xz ) -y x 

~V(B Z + ~f X y) 



V(B V + -f XZ ) 

V2lyz - B x 



1~V2 



V(B Z - ~f X y) 

V2lyz + B x 

-V _1_ ~ 1 ~'" 2 I -.' 1+V2 

Jxx 1 lyy 2 ' l zz 2 



(48) 



which goes over to Eq. (p3|) for v = 0, V2 = — 1. 
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The first-order average evolution operator, Eq. (26), 
contains the double integral of the commutator of the 
rotating-frame evolution operator T(U) taken at different 
moments U, i — 1,2. There are only a few combinations 
of trigonometric functions of the symmetrized angle ( [43] ) 
and the time-independent terms that result in non-zero 

contributions to T ^ for a single symmetric pulse with 
the duration r„: 



a = (sin((y9 — tp )}, 
a 2 = {sin(2tp-2tp')}, 
fx = (sin(2(^ - tp')), 



(49) 
(50) 
(51) 
(52) 
(53) 



where tp = tp(t), tp 1 
ages are defined as 



tp(t'), and the two-parameter aver- 



(f(t,t')) = - / dt / dt'f(t,t') 



(54) 



The parameters a and C completely characterize the Her- 
mitian evolution during the pulse; they were introduced 
and discussed in detail in Refs. 35, B2l 



B and decoherence matrix 7. 
pie cases, n = ±e x and n = 



We checked that in sim- 
±e„, the results coinside 



iden tically with those obtained directly as outlined in 
Sec. " 



IIIB1 



With these expressions, the evolution matrix Q for a 
sequence of pulses can be computed as a product of the 
evolution operators for individual pulses, 



Qi — Qo 



-p(O) 



p(l) 



2 1 ' 



G(r p 3 ) 



For a sequence of pulses, the net average evolution op- 
erator r soq can be obtained in terms of a logarithm of 
the obtained series (multiplied from the left by the cor- 
responding Qo,seq when it is not equal to 1). 

For the cycle 2s = XX of the Carr-Purcell sequence, 
the leading-order evolution operator is the same as for 
a single ^-pulse, Eq. fl48|), but with v — ► 0. The cor- 
responding expression in the subleading order is compli- 
cated and contains terms proportional to v, v 2 , vv-2, a, 
at2, C2; ah of these terms disappear in the limit of hard 
pulses. An antisymmetric version of the same sequence, 
2a = XX, results in the leading-order evolution opera- 



tor of the form (48), including the terms with v, but the 

subleading order disappears, r'' 1 (2a) = 0. Finally, for 
the phase-ramped four-pulse sequence 4a = XXXX, the 
leading-order effective decoherence operator is the same 
as for the Carr-Purcell cycle, 



2. Decoupling sequences of finite-length ix -pulses 

Given the computed terms of the average evolution op- 
erator for a single finite-duration pulse along the x axis, 
the corresponding expressions for a pulse along an ar- 
bitrary direction can be found by an orthogonal trans- 
formation, using the appropriately transformed vector 



r (0) (4a) = r (0) (2s) = r (0) (x) 



-(i) 



while the subleading order vanishes, Y (4a) = 0. 

For the two-dimensional decoupling sequence 4p = 
XYXY, the leading-order effective decoherence operator 
reads 



r (0) (4 P) 



/ lyy+lzz — 



ixx 4 



-S(B»+7»«) 



5f(S, + 7a «) 



Ixx + lz 



("~fxx 



l(lxy-B z ) 



iyy) 4 



(55) 



lz 



1+V2 

2 



Note that the corresponding diagonal terms are deter- 
mined by the diagonal terms of the original matrix 7 
and the parameter u 2 [see Eq. (^6|)], while off-diagonal 
terms are proportional to the effective width of the pulse 
v. For Hamiltonian evolution, 7 = 0, one could achieve 

r (4p) = by using self-refocusing pulses with v = 0. 
The next order term is complicated, but it is eliminated 
for the antisymmetric sequence, 8a = XYXYYXYX. 
Specifically, for this sequence we obtain the same leading- 
order expression r (0) (8a) = r (0) (4p) [Eq. (fBI)], and 



r 1 (8a) = 0. On the other hand, for the symmetric 
sequence 8s = XYXYYXYX, it is the leading-order off- 
diagonal corrections proportional to v that vanish, 



r (0) (8s) 



r (0) (4 P)! 



while the subleading order contains terms proportional 
to v, a, ot2, and £2. All of the corrections can be elimi- 
nated by a supercycle of phase ramping which leads to an 
antisymmetric 16-pulse sequence 16a constructed as the 
8s sequence followed by the same sequence with inverted 
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pulse 


00 


V 


V2 


a/2 


Q2/2 


c 


C2 


A 4 


<M(t - 




7>/2) 




00 


cos ^ 


COS 00 


1 • 1 
g sin 0o 


1 ■ nl 

g sin 20o 


I sin 3a 


1 • i 
4 sin 0o 


1 ■ 3d>n 


n5{t- 


t p 


/2) 




7T 





-1 








1/4 





-1/4 


Go.oi [ 


5 


J ] 






7T 


0.0211 


-0.9709 


0.0104 


0.000047 


0.24996 


0.00023 


-0.2354 


Go. 10 


5 


il 






7T 


0.2107 


-0.7086 


0.0872 


0.0047 


0.2458 


0.0233 


-0.1035 


Si 


19 








7T 





-0.6135 


0.0333 


0.0393 


0.2382 


-0.0737 


-0.1020 


S 2 


19 








7T 





-0.6675 


0.0250 


0.0298 


0.2414 


-0.0557853 


-0.1171 


Qi 


1£ 








7T 





-0.6761 





-0.0079 


0.2399 


0.0027 


-0.1234 


Q 2 


1£ 








7T 





-0.7138 





-0.0065 


0.2422 


0.0022 


-0.1342 


Wii 


( 


7T 


) 






7T 








0.0511 


-0.0039 


0.1884 


-0.1014 


0.0353 


Wia( 


7T 


1 






7T 








0.0400 


-0.0164 


0.1904 


-0.0871 


0.0413 


Wai( 


7T 


) 






7T 











0.0088 


0.0072 


0.0677 


-0.0093 


W 22 


( 


7T 


) 






7T 











0.0107 


0.0634 


0.0415 


-0.0035 


W 3 i( 


7T 


) 






7T 











0.00061 


0.0436 





0.0014 


W 32 ( 


7T 


) 






7T 











0.00046 


0.0847 





0.0146 


Fj 










7T 


0.0018 


0.3307 


0.0237 


-0.01018 


0.1134 


-0.0260 


0.0680 


PKRUi( 


7T 


)[ 






7T 





-0.6880 


0.0278 


0.0368 


0.2420 


-0.0614 


-0.1250 


PKRU 2 ( 


7T 


)[ 






7T 





-0.1501 





-0.0078 





0.0866 


-0.0047 


f5(i- 


Tp/2) 




tt/2 


V2/2 





1/8 





\/2/8 


1/4 


\/2/8 


Go. 01 


(77/2) 




tt/2 


0.7136 


0.0211 


0.1272 


0.0104 


0.1767 


0.2500 


0.1872 


Go. 10 


(tt/2) 




tt/2 


0.7722 


0.2107 


0.1388 


0.0872 


0.1706 


0.2458 


0.2599 


Qi(tt/2) 


52 


j 




tt/2 





-0.2906 





-0.0023 


0.2021 


0.0040 


-0.0508 


<?2(7r/2) 


52 


j 




tt/2 





-0.0109 





0.0035 


0.1617 


0.0730 


0.00077 


Wu( 


7T 


/2) 






tt/2 








0.0106 


-0.0022 


0.1787 


0.0114 


0.0193 


Wia( 


7T 


/2) 






tt/2 








-0.0057 


-0.0020 


0.1756 


0.0482 


0.0103 


Wai( 


7T 


/2) 






tt/2 











-0.0059 


0.1796 


0.0301 


0.0190 


W 22 ( 


7T 


/2) 






tt/2 











-0.0021 


0.1771 


0.0324 


0.0129 


PKRUi( 


TT/2) 


54 


I 


tt/2 





-0.2529 


-0.0260 


-0.0242 


0.1992 


0.0541 


-0.0422 


PKRU 2 ( 


7T 


/2) 


54 


I 


tt/2 





0.2711 





0.0127 





0.0483 


0.00028 



TABLE I: Expansion coefficients for different pulse shapes. See text for the definitions. 



pulses. Specifically, 



inverse decoherence time becomes 



r (0) (i6a) = r (0) (8s) = r (0) (8a) 

r (1) (16a) = f (1) (8a) =0, 



(56) 
(57) 



where r (0) (8a) 



-(0), 



r v "'(4p) is given by Eq. @. Wo 
should mention that for Hamiltonian evolution, 7 = 0, 



achieve r (0) (8s) 




one could 

second-order pulses with v — a = |35|, 15 

Comparing these results with the analogous sequences 
of hard pulses, we see that an equivalent cancellation 
with soft pulses requires careful pulse shaping or dou- 
bling the number of pulses in the sequence. On the other 
hand, while sequences of hard tt pulses do not modify the 
diagonal part of the effective decoherence operator, this 
is not so with shaped pulses where 112 / 0. This corre- 
sponds to a reduction of the effective decoherence rate 
T 1 due to redistribution between the channels. For ex- 
ample, with the NMR decoherence model, Eq. (|i~7j), the 



T l 1 = 2 7- {l ~ 14>) 1 



(58) 



for 7 > 70 and limited control fields (\v%\ < 1), this 
is smaller than the original T^ 1 = 27. Note that the 
reduction is achieved solely by redistribution of decoher- 
ence between the channels, that is, by a corresponding 
increase in T 2 ■ For full symmetrization [Eq. (f28|)l one 



using needs v-i = 1/3. 



3. Sequences of finite-length tt/2 pulses 

For the family of sequences 12, 24, 48 [sec Eqs. ( |35| ) - 
([37])] of 7r/2 pulses, the results are similar to those with 
hard pulses. Specifically, when used with NMR decoher- 
ence model (|l7|), the sequence 12 achieves leading-order 
coherence matrix symmetrization but not the decoupling. 
The symmetrized (and phase ramped) sequences 24 and 
48 achieve symmetrization and decoupling in the leading 
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and subleading orders, respectively. More explicitly, for 
symmetric pulse shapes and the decoherence model (17), 



r (0) (24) = r (0) (48) = |(27+ 7 *)1, T (1) (48) = 0. (59) 



The non-zero matrix elements in T 1 (24) contain terms 
scaling as the products of (7 — 70) and the external field 
components in all combinations; all of the coefficients 
cannot be eliminated merely by pulse shaping. In the 

special case = 0, the correction I^ 1 ' > (24) oc (7 — 7^) 2 
is suppressed for pulse shapes with V2 = 0L2 = 0. 

With generic decoherence matrix 7, r^ 1 (48) = for 
90-degree pulses with v = V2 = 0. In this case the 
full symmetrization of the decoherence matrix can be 
achieved as long as ^ xx = ^ vv , which can be always made 
to be the case by an appropriate choice of the basis in 
the x-y plane. 



C. Pulse shaping 

The obtained analytical results imply two possible ap- 
plications for pulse shaping. First, a pulse can be shaped 
to have the coefficients v, a, a2, and C2 zero, so that the 
second-order correction to effective decoherence operator 
is zero already for the 8-pulse sequence 8s. The resulting 
effective decoherence operator T will be diagonal, with 
the matrix elements determined by the diagonal terms of 
the original decoherence matrix 7 and the parameter w 2 
[cf. the diagonal elements in Eq. (55)]. We note that up 
to terms of second order in pulse duration, one might as 
well use the longer sequence 16a to cancel the terms in 
r proportional to coefficients v, a, a-z, £2- 

The value V2 = — 1 corresponds to an ideal hard inver- 
sion pulse, in which case the diagonal matrix elements of 
the matrix T are not modified. On the other hand, for 
V2 = 1/3, assuming j xx = ^yy (which can always made to 
be the case by a rotation in the x-y plane), the effective 
decoherence operator is fully symetrized [see Eq. (p8|)]. 
This is the second potential application for pulse shap- 
ing. Again, to second order, the same symmetrization 
can be also achieved with the composite-pulse sequence 
48 using 7r/2-pulses with v — V2 = (generic decoher- 
ence matrix) , or any pulses for the special case of NMR 
decoherence model (17). 

We followed Refs. IE , ^ to construct a number of new 
it/2 and inversion (0o = f) pulse shapes with v — V2 = 
[Wi s (<f> )], v = v 2 = a = [W 2lI (0o)], and v = v 2 = 
a = (2 = [W^s^o)], where the parameter s = 1,2 
determines the number of derivatives that vanish at the 
ends of the pulse interval (2 and 4, respectively). The 
parameters of the new and previously constructed pulses 
are listed in Tab. [j], and their coefficients in Tab. O. The 
shapes of the pulses we used in the simulations are shown 
in Fig. §. 




FIG. 2: (Color online) Pulse shapes used in the simulations. 
The Gaussian shape G x (y) = [(27r) 1/2 2;r p ]" 1 exp(-0.5 y 2 /x 2 ) 
with x = 0.10 and y = t/r p — 0.5 is \/2 wider than that used 
in Ref. IB 511 . The shapes Si, Qi were originally constructed 
in Ref. |19f . The shapes F\, Wn, W21 are constructed in this 
work, see Tab. [n| for the coefficients. 



IV. SIMULATIONS 

Our analytical analysis was only done to second or- 
der cumulant expansion, and assuming time- independent 
fields The results are asymptotically true in the 

limit where the sequences are sufficiently short, that is, 
tB C 1, 7T C 1, t C r c , where r c ~ uj^ 1 is the cor- 
relation time of the slow part of the bath. To check for 
possible effect of ignored terms, we performed numeri- 
cal simulations of the controlled dynamics of a qubit in 
the presence of time-dependent correlated classical Gaus- 
sian fields B^ = Bfj,(t), as well as Markovian decoher- 
ence described by the Lindblad equation (11), (12). We 



specifically used the NMR decoherence model (17) in the 
regime dominated by phase fluctuations, 7^ > 0, 7 = 
[note that with a non-zero 7 the decoherence matrix can 
be decomposed as 7 = 7I + (7^ — 7) diag(0, 0, 1); the 
part proportional to the identity matrix commutes with 
decoupling pulses and can be trivially eliminated.] The 
fields B/j.(t) were chosen as zero- average correlated Gaus- 
sian random functions with the correlation function 



(B^(t)B v (t')) = 5^B$g(t-t'), g(t) = e 



-t 2 /2-r c 2 



(60) 



and the fixed value of the correlation time r c = 8r p . The 
control fields V^(t), [i = x,y, were generated according to 
the chosen decoupling sequence and the pulse shape, with 
the fc-th pulse of the sequence fitting in the interval (k — 
1)t p < t < kr p ; the corresponding decoupling period r = 
nTp is given by the number n of pulses in the sequence. 
For every sequence, pulse shape, and the realization of 



the random fields B p (t), we solved Eqs. (p.5[), (|16| ) with 
TZ = for three sets of initial conditions with i? M = 1, 
fj, = x,y, z; the solutions correspond to the columns of the 
evolution matrix Q(t) [Eq. (|2l])]. The average decoupling 
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pulse 


AO 


Al 


A2 


A3 


A4 


A5 


A6 


A7 


Fi 


0.5 


-1.419474 


-2.048028 


1.549555 


1.435813 


-0.017867 








0.5 


-1.242022 


-1.009075 


0.700828 


0.530624 


1652161644 


0.277982 


0.241663 


Wi2(n) 


0.5 


-1.291342 


-0.753726 


1.499438 


0.364546 


0.012680 


-0.069983 


-0.261614 


W 2 l(7T) 


0.5 


3.056086 


-1.295369 


-1.689687 


-0.062202 


-0.366646 


-0.142183 




W22(k) 


0.5 


2.776007 


-2.473314 


-1.782314 


0.958211 


-0.444991 


0.300165 


0.166236 


W 3 l(7r) 


0.5 


-1.110710 


-3.692547 


1.248118 


0.990698 


1.394824 


0.669618 




W32(n) 


0.5 


-1.686664 


-2.108402 


3.362253 


1.029286 


-0.260405 


-0.836068 




Wuin/2) 


0.25 


2.011311 


0.041292 


1.381531 


0.262448 


0.076040 






W 12 (ir/2) 


0.25 


2.023581 


0.920572 


1.341484 


-0.113434 


-0.144034 


-0.231008 




W 21 (it/2) 


0.25 


2.018463 


0.588295 


1.393403 


-0.206226 


0.095943 


-0.1029524 




W 2 2(it/2) 


0.25 


2.018283 


0.608538 


1.386685 


0.088935 


0.024615 


-0.134584 


-0.205904 



TABLE II: Pulse coefficients for various shapes constructed for this work. The control field during the pulse, < t < t p , is 
represented as V(t) = 2tt A n cos(27rnt/r p ). 



fidelity ( J29| ) was then evaluated using the expression 

W,» = i(l + ^>). (61) 

at the time moments t s commensurate with the sequence 
duration t s = st — snT p . 

The calculation results for a particular sample of the 
fields Bfj,(t) and several sequences of inversion (<jf>o = it) 
Gaussian pulses are shown in Fig. ||. The line marked as 
"ideal" corresponds to the best achievable fidelity 

(F(t)) idcal = i(2 + e^ t ); (62) 

with the chosen value 7^ = 27r x 10~ 3 / r p at the end of the 
simulation we get (F(512T p ))idoai ~ 0.680. The points 
marked "no control" corresponds to average fidelity in 
the absence of decoupling pulses; the coherence is lost 
after just a few pulse durations. While decoupling is not 
particularly efficient (and noise is evident) for the 4-pulse 
sequence (set of points marked 4p in Fig. ||) , the fidelity 
is improved and noise markedly reduced for the sequence 
8s, and even more so for the sequence 16a. In fact, the 
points for sequence 16a are very close to the dashed line 
which represents the leading-order contribution at B = 0, 
see Eq. ( J63| ) below. 

In the presence of the Markovian decoherence, we can 
view the average infidelity as composed of three terms. 
First is the inavoidable infidelity due to the Markovian 
decoherence alone, in the absence of any control fields or 
low-frequency noise. For the case of pure dephasing, it 
is given by 1 — (F(t))idcai, see Eq. (|62]). For the value of 
7^ used in our simulation, (.F(i)}ideai is shown in Fig. || 
with the solid line. Second is the fidelity loss due to the 
redistribution of the original decoherence rate(s) over the 
channels. In our simulations, for the sequences 4p, 8s, 
16a of 7r-pulses other than Go. 10, this part of the average 
infidelity is to a very good accuracy determined by the 
coefficient v 2 of the corresponding pulses, 

{F(t)) = \ (3 + e-^»«'+2e-T 2 »* t ), (63) 




0.3 1 ' 1 ' 1 !-l 

100 200 300 400 500 



FIG. 3: (Color online) Single-qubit average fidelity for a par- 
ticular realization of time-dependent fields and Markovian de- 
phasing with the rate 7^ = 2ty x 10 _3 /r p . Open boxes con- 
nected with the dashed line marked "no control" correspond 
to the absence of decoupling. Other symbols are numerical 
data for the indicated decoupling sequences of Gaussian inver- 
sion pulses. The black solid line marked "ideal" corresponds 
to maximum achievable average fidelity. Dashed line corre- 
sponds to Eq. ( |63| ) with v 2 = —0.7086 as appropriate for the 
pulse Go. 10 [see Tab. [j]. 

7icff = 14, — g — ' 72cfi = 70 — ^ — . (64) 

see Tab. | and Fig. |. 

Finally, the third part of the average infidelity is en- 
tirely due to the presence of the low-frequency random 
fields; it is defined as the difference between the aver- 
age fidelity of the controlled system in the absence of 
low-frequency fields, and that in their presence, AF(t) = 
(F(t))o — (F(t)). It is this quantity that directly charac- 
terizes the ineffectiveness of the decoupling against the 
low- frequency fields. 

The corresponding plots are shown in Figs. || (sequence 
8s) and |^ (sequence 16a), as well as in Figs. (?], || with 



11 




FIG. 4: (Color online) Single-qubit average fidelity loss due 
to redistribution of the decoherence rates over directions in 
the absence of low-frequency fields. Symbols are numerical 
data for pulses and sequences as indicated, and lines are com- 
puted analytically using Eq. ( |63| with corresponding V2 from 
Tab. [l| This contribution to infidelity would be zero for hard 
(S- function) pulses with V2 = — 1. 



FIG. 5: (Color online) Single-qubit decoupling error for the 8- 
pulse sequence 8s with several pulse shapes, as indicated. The 
decoupling error is small for the second-order pulse W / 2i(7r) 
and, somewhat surprisingly, for the symmetrizing shape F%. 
With 7^, = 0, the infidelity for this pulse shape with sequence 
8s grows similarly to that of first-order pulses Si and Wu 
[Fig.0. 



70 = for comparison purposes. The general trend is 
consistent with the expectations based on analytical ex- 
pansion. For example, with the Gaussian pulses we get 
the largest error, with or without the Markovian dephas- 
ing; the decoupling error is smaller with the 16a sequence 
[as compared with the 8s] where the corrections including 
the subleading order are suppressed. Similarly, specially 
designed pulse shapes with a larger number of suppressed 
parameters lead to improved performance in the presence 
of Markovian dephasing: compare, e.g., the curves for 
pulses Si, Wu(tt), and W^iW- 

One marked exception is the pulse shape F\ , for which 
v-2 ~ +1/3 while the other coefficients have magnitude 
comparable to those for, e.g., the shape Si (apart from 
the coefficient v which is zero for the shape Si). While 
in the absence of Markovian dephasing, the decoupling 
error for the pulse F\ is qualitatively similar to that for 
the lst-order shapes [Figs. ^, with 7^ ^ the decou- 
pling error is actually closer to that for the second-order 
pulses with the sequence 8s [Fig. || and is the smallest of 
the shapes we considered for the sequence 16a [Fig. §. 
We believe that this, at least in part, is associated with 
the symmetrization of the effective decoherence tensor 
achieved in the presence of the pulse Fi , which renders 
the resulting Markovian decoherence (effectively, the am- 
plitude damping) independent from the decoupling, and 
thus reduces the contribution from high orders not in- 
cluded in our analytical calculations. 

In comparison with the sequences of 7r-pulses, the ba- 
sic unit in the sequences of n/2 pulses, 12p, 24s, and 
48s, is three times longer. As a result, weaker random 
fields B^(t) with longer correlation times are required 
to achieve comparable decoupling accuracy. In addition, 
with longer coherence time, longer evolution times are 



required to separate the effects of the transients near the 
beginning of the simulation interval. By these reasons, 
we do not discuss the performance of these sequences in 
detail. For the parameters as in Figs. |^, [| the decoupling 
error with the sequences 24s and 48s saturates at around 
twenty times that for the sequences 8s, 16a, while with 
t c = 32t p and the same Bq = 0.1/t^, the decoupling er- 
ror is comparable to that in Figs, ^. Note that these 
decoupling errors are still small compared to the fidelity 
loss due to redistribution of the decoherence rates over 
directions, which for these sequences is identical to that 
of sequences 8s, 16a with the symmetrizing pulse Fi [see 

Fig. a. 



V. CONCLUSIONS 

We considered the dynamical decoupling in a simple 
decoherence model simulating the presence of both low- 
and high-frequency environment modes [Fig. |], having 
in mind applications of DD in combination with QECCs. 
We modeled the effect of low-frequency degrees of free- 
dom in terms of classical correlated noise, and that of 
the fast degrees of freedom with the help of the Marko- 
vian master equation in the Lindblad form, see Eqs. ji"l]), 
(1C). The combined effect is the non-Hermitian evolution 
(15) of the Bloch vector with the instantaneous decoher- 
ence operator ( |l6| ) [see Eq. ( |l7| ) for the special case of 
NMR decoherence model]. In the presence of dynam- 
ical decoupling with the sequence period r, the values 
of the average qubit fidelity at the commensurate time 
moments t s = st are determined by the effective deco- 
herence matrix T [Eq. (|2^)], an analogue of the average 
Hamiltonian. The trace of the effective decoherence ma- 
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FIG. 6: (Color online) As in Fig. ^ but for 16-pulse sequence 
16a. For this sequence, and with time-independent fields _B M , 
there would be no contribution from the first two orders of 
the effective decoherence operator [Eqs. jjil), (p7|)1; the decou- 
pling error is determined by the non-adiabaticity and by the 
terms of higher order. Generic pulse shapes Go.io, Si, and Qi 
show the largest decoupling error, while it is the smallest for 
the symmetrizing shape Fi (technically, zeroth order, 
V2 ~ 1/3). For shapes F\, Wn, and W21, the decoupling error 
here is actually smaller than the infidelity in Fig. pi We be- 
lieve this is an artefact of the definition of the decoupling error 
related to the full average fidelity approaching saturation at 
1/2 [Fig. § . 



trix is conserved, see Eq. (|2j); in agreement with general 
expectations, this implies that dynamical decoupling can 
only decrease the fidelity in the presence of Markovian 
degrees of freedom alone. 

This is always the case for DD with soft pulses, which 
necessarily leads to redistribution of the decoherence 
rates over the directions, leading to some fidelity reduc- 
tion even in the absence of the slow degrees of freedom. 
For evolution time small compared to the decoherence 
time ~ 7 _1 (as required for efficient error correction), 
this reduction is a relatively small effect. On the flip side, 
this redistribution causes symmetrization of decoherence 
operator which reduces the effect of the non-Hamiltonian 
dynamics associated with the fast degrees of freedom on 
the decoupling accuracy. For example, if the symmetriza- 
tion ( plf ) is achieved at the end of a basic decoupling 
cycle, decoherence is expected to have no effect on addi- 
tional cancellations achieved in a supercycle obtained by, 
e.g., phase-ramping the basic cycle. 

To analyze relative importance of these effects, we con- 
sidered several decoupling sequences of both hard and 
soft (generic 1st- and 2nd-order and specially shaped) 
pulses which lead to various degrees of symmetrization 
of the decoherence operator. For such sequences, in the 
static limit where the slow degrees of freedom become 
time-independent, we constructed analytically the first 
two terms of the (cumulant) expansion of the effective 
decoherence operator T in powers of the sequence pe- 




FIG. 7: (Color online) As in Fig. |s| but in the absence of 
Markovian dephasing, 7^ = 0. Infidelity curves with 2nd- 
order pulses Qi and W21 have a finite slope, unlike the plots 
in Ref. ^ where a smaller r.m.s. value Bo of the slow fields 
was used. We checked that the slopes for the 2nd-order pulse 
shapes are greately reduced if Bo is reduced by a factor of two 
(not shown). 
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FIG. 8: (Color online) As in Fig. f?| but for the sequence 16a. 
The slopes for all the pulse shapes excluding Go.io are greately 
reduced if Bo is reduced by a factor of two (not shown). For 
the Gaussian pulses Go. in, the second-order cancellation is 
achieved only after the complete sequence with the duration 
r = 16r p ; the slope reduction happens for slower fields B M (t) 
with t c — 16r p (not shown). 

riod r. With the help of pulse shaping and/or phase 
ramping, we could ensure that the leading-order average 

decoherence operator ' is diagonal and independent 
of the components of the slow field, while the subleading 

contribution disappears identically, r 1 = 0. Results of 
Refs. |lj| ^l] suggest that such sequences might also be ef- 
fective in the presence of sufficiently weak slowly- varying 
fields. 

Numerical simulations generally confirmed these ex- 
pectations. In the studied parameter range, the best 
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decoupling accuracy is achieved with the sixteen-pulse 
sequence 16a. With time-independent external fields 
the sequence provides complete decoupling of the slow 
fields in the first two orders of the cumulant expansion 
independent of the shape of symmetric inversion pulses, 
see Eqs. ([So]), (|57|). The effects of higher order terms 
and non-adiabaticity of the fields B^{t) are captured in 
the simulations, see Sec. IV. The best performance is 
achieved with the "symmetrizing" pulse shape F\ con- 
structed with the only requirement that the diagonal part 
of the leading-order effective decoherence matrix (|5^) for 
the 4-pulse sequence 4p is symmetric, i.e., vi ps +1/3, 
see Sec. IIIB 1 and Tab. In fact, in our simulations the 
part of the coherence loss associated with the slow fields 
is smaller than that in the absence of the high-frequency 
degrees of freedom, cf. Figs. ||, || and 0, ||. Overall, the 
suppression of the decoherence due to slow fields associ- 
ated with the dynamical decoupling increases the coher- 



ence time by orders of magnitude. 

We thus demonstrated in principle the effectiveness of 
DD in suppressing the effect of low-frequency degrees of 
freedom in the presence of high-frequency modes which 
cannot be eliminated by dynamical decoupling. This 
opens up possible applications in combined coherence 
protection techniques concatenating dynamical decou- 
pling with QECC at higher levels. We postpone the cor- 
responding discussion to a futhcr publication fl56| . 



Acknowledgments 

We are grateful to Sasha Korotkov, Daniel Lidar, Len 
Mueller, and Gotz Uhrig for useful discussions. This 
research was supported in part by the NSF grant No. 
0622242. 



[1] C. P. Slichter, Principles of Magnetic Resonance, 
Springer- Verlag, New York, 3rd edition, 1992. 

[2] P. Hodgkinson and L. Emsley, Prog. NMR Spectr. 36, 
201 (2000). 

[3] L. M. K. Vandersypen and I. L. Chuang, Reviews of Mod- 
ern Physics 76, 1037 (2004). 

[4] U. Haeberlen and J. S. Waugh, Phys. Rev. 175, 453 
(1968). 

[5] U. Haeberlen, J. S. Waugh, Phys. Rev. 185, 420 (1969). 
[6] A. Pines, J. S. Waugh, J. Mag. Res. 8, 354 (1972). 
[7] A. Pines, J. S. Waugh, Phys. Lett. A 47, 337 (1974). 
[8] L. Viola and S. Lloyd, Phys. Rev. A 58, 2733 (1998). 
[9] D. Vitali and P. Tombesi, Phys. Rev. A 59, 4178 (1999). 
[10] L. Viola, E. Knill, and S. Lloyd, Phys. Rev. Lett. 82, 
2417 (1999). 

[11] L. Viola, S. Lloyd, and E. Knill, Phys. Rev. Lett. 83, 
4888 (1999). 

[12] L. Viola, E. Knill, and S. Lloyd, Phys. Rev. Lett. 85, 
3520 (2000). 

[13] D. Vitali and P. Tombesi, Phys. Rev. A 65, 012305 
(2002). 

[14] C. Uchiyama and M. Aihara, Phys. Rev. A 66, 032313 

(2002) . 

[15] M. S. Byrd and D. A. Lidar, Phys. Rev. A 67, 012324 

(2003) . 

[16] K. Shiokawa and D. A. Lidar, Phys. Rev. A 69, 

030302(R) (2004). 
[17] L. Faoro and L. Viola, Phys. Rev. Lett. 92, 117905 

(2004) . 

[18] P. Facchi et al., Phys. Rev. A 71, 022302 (2005). 
[19] P. Sengupta and L. P. Pryadko, Phys. Rev. Lett. 95, 
037202 (2005). 

[20] K. Khodjasteh and D. A. Lidar, Phys. Rev. Lett. 95, 
180501 (2005). 

[21] L. P. Pryadko and P. Sengupta, Phys. Rev. B 73, 085321 
(2006). 

[22] P. Chen, Phys. Rev. A 75, 062301 (2007). 
[23] G. S. Uhrig, Phys. Rev. Lett. 98, 100504 (2007). 
[24] T. E. Hodgson, L. Viola, and I. DAmico, Phys. Rev. B 
78, 165311 (2008). 



[25] S. Pasini, T. Fischer, P. Karbach, and G. S. Uhrig, Phys. 

Rev. A 77, 032315 (2008). 
[26] M. S. Byrd and D. A. Lidar, J. of Modern Opt. 50, 1285 

(2003). 

[27] L. Viola, Phys. Rev. A 66, 012307 (2002). 

[28] M. Mottonen, R. de Sousa, J. Zhang, and K. B. Whaley, 

Phys. Rev. A 73, 022332 (2006). 
[29] P. Kuopanportti et al., Phys. Rev. A 77, 032334 (2008). 
[30] L. Cywihski, R. M. Lutchyn, C. P. Nave, and S. D. 

Sarma, Phys. Rev. B 77, 174509 (2008). 
[31] P. Rebentrost, I. Serban, T. Schulte-Herbriiggen, and 

F. K. Wilhelm, Phys. Rev. Lett. 102, 090401 (2009). 
[32] G. S. Uhrig, New J. Phys. 10, 083024 (2008). 
[33] J. S. Waugh, L. M. Huber, and U. Haeberlen, Phys. Rev. 

Lett. 20, 180 (1968). 
[34] J. S. Waugh, C. H. Wang, L. M. Huber, and R. L. Void, 

J. Chem. Phys. 48, 652 (1968). 
[35] L. P. Pryadko and G. Quiroz, Phys. Rev. A 77, 012330/1 

(2007) . 

[36] O. V. Konstantinov and V. I. Perel, Zh. Exp. Teor. Fiz. 

39, 197 (1960). 
[37] E. B. Davies, Comm. Math. Phys. 39, 91 (1974). 
[38] M. I. Dykman, Physica Status Solidi (b) 88, 463 (1978). 
[39] R. Alicki, Phys. Rev. A 40, 4077 (1989). 
[40] A. G. Kofman and G. Kurizki, Phys. Rev. Lett. 93, 

130406 (2004). 
[41] G. Lindblad, Commun. Math. Phys. (1976). 
[42] R. Alicki, D. A. Lidar, and P. Zanardi, Phys. Rev. A 73, 

052311 (2006). 

[43] I. Gourlay and J. F. Snowdon, Phys. Rev. A 62, 022308 
(2000). 

[44] L. Ioffe and M. Mezard, Phys. Rev. A 75, 032345 (2007). 
[45] Z. W. E. Evans, A. M. Stephens, J. H. Cole, 

and L. C. L. Hollenberg, unpublished, preprint 

arXiv:0709.3875 (2007). 
[46] P. K. Sarvepalli, A. Klappenecker, and M. Rotteler, Proc. 

R. Soc. A 465, 1645 (2009). 
[47] P. Aliferis and J. Preskill, Phys. Rev. A 78, 052331 

(2008) . 

[48] A. M. Stephens, Z. W. E. Evans, S. J. Devitt, , and 



14 



L. C. L. Hollenberg, Phys. Rev. A 77, 062335 (2008). 
[49] P. Zanardi, Phys. Lett. A 258, 77 (1999). 
[50] W.-K. Rhim, D. D. Elleman, and R. W. Vaughan, J. 

Chem. Phys. 59, 3740 (1973). 
[51] G. A. Morris and A. Gibbs, Magnetic Resonance in 

Chemistry 29, 83 (1991). 
[52] L. P. Pryadko and P. Sengupta, Phys. Rev. A 78, 032336 

(2008). 



[53] W. S. Warren, J. Chem. Phys. 81, 5437 (1984). 

[54] S. Pasini, P. Karbach, C. Raas, and G. S. Uhrig, unpub- 
lished, preprint arXiv:0906.2070 (2009). 

[55] C. Bauer, R. Freeman, T. Frenkiel, J. Keeler, and A. J. 
Shaka, J. Mag. Res. 58, 442 (1984). 

[56] Y. Li, D. Lidar, and L. P. Pryadko, In preparation. 



